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Abstract 

66" 
On ■ 
G\ : 

We develop a first principles, microscopic theory of impurity atom scat- 
\ tering from inhomogeneous quantum liquids such as adsorbed films, slabs, or 

Ph ' clusters of 4 He. The theory is built upon a quantitative, microscopic descrip- 

tion of the ground state of both the host liquid as well as the impurity atom. 
Dynamic effects are treated by allowing all ground-state correlation functions 
t— 1 , to be time-dependent. 

^ I Our description includes both the elastic and inelastic coupling of impu- 

£Sj ' rity motion to the excitations of the host liquid. As a specific example, we 

study the scattering of 3 He atoms from adsorbed 4 He films. We examine the 
q ■ dependence of "quantum reflection" on the substrate, and the consequences of 

impurity bound states, resonances, and background excitations for scattering 
^ I properties. 

A thorough analysis of the theoretical approach and the physical circum- 
2 stances point towards the essential role played by inelastic processes which 

■ determine almost exclusively the reflection probabilities. The coupling to im- 

purity resonances within the film leads to a visible dependence of the reflection 
q I coefficient on the direction of the impinging particle. 

> : 

• i-H , 

X ' I- INTRODUCTION 



Dynamic scattering processes of helium atoms from low temperature liquid 4 He films and 
the bulk fluid in the vicinity of a free surface continue to be a subject of considerable interest. 
Experimental information is available mostly for 4 He scattering processes, connected with 
quantum reflection and quantum evaporationiH 4 !, as well as the surface reflect ivityiHi'0. Due 
to experimental difficulties, there are o„lv few data for *He scattering but there is also 
interest (experimental0Hll, and theoretica0~@ ) in the dynamics of atomic Hydrogen atoms 
on 4 He surfaces for which our theory also applies. 

This paper follows up on a line of work studying the properties and the dynamic fea- 
tures of quantum liquid films from a manifestly microscopic point of view. Most relevant 
for the present work are papers designing the theory for the background host liquidE, its 
excitations0'@, and the dynamics of atomic impurities^. In that work, we have used the 
method of correlated variational wave functions which has in many situation proven to be a 
computationally efficient, precise, and robust method for the purpose of studying strongly 
interacting quantum liquids. Even the simplest approximation of the theory has in the 
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past eiven quite satisfactory results on the nature of the impurity states!!, their effective 
massE3 and the impurity-impurity interaction!! in inhomogeneous geometries. The reason 
for the qualitative success of the theory is that it contains a consistent treatment of both 
the short- and the long-range structure of the system. This implies that both the low- and 
the high-lying excitations are treated accurately. 

The present paper complements a similar study of the scattering of 4 He atoms from 4 He 
slabs@; the problem at hand is somewhat simpler since there is no need to fully symmetrize 
the wave function of the background system and the impinging particle. Another major 
physical difference to the scattering of 4 He particles is that in the latter case one might 
observe@~i3 the coupling to the Bose-Einstein condensate, whereas in the present case one 
can couple both to phonon-like and to single particle excitations. Nevertheless we will see 
that many similarities exits between the two problems: The scattering process is dominated 
by inelastic channels, mostly the coupling to ripplonic excitations. 

Generally, the impinging particle can, in the presence of other particles like the film of 
4 He under consideration here, scatter into three types of channels: 

1. Elastic reflection: The incoming particle, characterized by the wave vector (kiK&jJ, 
is elastically reflected with a probability |-R| 2 . It creates virtual excitations of the 
background, but transfers no energy. 

2. Inelastic scattering: with a probability r ine i the particle loses some energy to an exci- 
tation of the film, and retains enough energy to leave the attractive potential of the 
film and the substrate. The film excitation can be either a collective wave (ripplon, 
phonon), or a single 4 He that is elevated above the chemical potential /14 and leaves the 
film. The creation of several excitations is in principle also included in our theoretical 
description, but it is ignored in the linearized treatment of the equations of motion. 

3. Adsorption: as in the previous case, the film is excited, but the particle is adsorbed to 
the film. The corresponding sticking coefficient s is the probability for this process. 

These three types of processes are depicted in Fig. [l]. Because of the hermiticity of the 
many-body Hamiltonian for N 4 He atoms and the 3 He impurity, we have 

\R\ 2 + r ine i + s = 1. (1.1) 

This work focusses on the calculation of elastic scattering because the impinging particle 
couples, in particular at low energies, predominantly to the low-lying, bound excitations of 
the background film and the impurity atom. We shall argue below that, basically for phase- 
space reasons, inelastic processes are expected to be less important than either elastic, or 
total absorption processes. 

Since most of the theoretical tools of the present study have been derived in Ref. |20], we 
outline in Sec. |I| only briefly the theoretical methods and the basic equations to be solved. 
The scattering problem will be be formulated in terms of a non-local, energy dependent 
"optical potential" which depends explicitly on the coupling of the impinging particle to 
background and impurity excitations. 

The results of our calculations are discussed in section |V|. To cover a variety of physical 
situations, we will present results for several of the systems that were studied extensively in 
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our previous calculations: These will range from strongly bound films on a model graphite 
substrate that is covered with two layers of solid helium, to a very weakly bound model, 
described by a rather thick, metastable film on a Cesium substrate. We first discuss the 
possible excitations of the background systems, and then present results for the surface 
reflectivity as a function of impact energy and angle for some of those systems. At very low 
energies, we will encounter the effect of "quantum reflection"ilHi3IllEBEil; with increasing 
impact energies we also can analyze the influence of surface excitations (ripplons) and the 
Andreev state, phonon/roton creation, and under certain circumstance the coupling to an 
"Andreev resonance" of the impurity particle close to the substrate. 



II. MICROSCOPIC THEORY 

The theoretical description of 4 He films and impurity properties starts with a description 
of the ground state of the background system. Next, a single impurity is added, and finally 
this impurity is allowed to move. The technical derivation and in particular the important 
verification of our theoretical tools have been presented in a series of previous papers§10@, 
we will therefore discuss the theoretical background only briefly. 



A. The Background Liquid 



In the first step, one calculates the properties of the background helium film. The only 
phenomenological input to the theory is the microscopic Hamiltonian 



Hs= E 

Ki<N 



h 2 



2m 



B 



+ E v-(ki-'il), 

l<i<j<N 



(2.1) 



where V(|r^ — r^l) is the 4 He- 4 He interaction, and Z7 S ub( r ) is the external "substrate" potential. 
The many-body wave function is modeled by the Jastrow-Feenberg ansatz 



*j\r(ri, . . . ,rjy) = exp 



E M r i)+ E «a( r *> r i) + E u 3 (ri,rj,r k ) 

l<i<N l<i<j<N l<i<j<k<N 



• (2-2) 



An essential part of the method is the optimization of the many-body correlations by solving 
the Euler equations 



5E 



N 



,ri,...,r r , 







[n = 1, 2, 3) 



(2.3) 



where En is the energy expectation value of the A-particle Hamiltonian (|2.1|) with respect 
to the wave function 



E 



N 



J d?ri . . . d 3 r N ^ n (t u t^Hn^n^, . . . , r N ) 
f d 3 ^ . . . d 3 r N ty 2 N (r u . . . ,r N ) 



(2.4) 



The energy is evaluated using the hyper-netted chain (HNC) hierarchy of integral 
equations^; "elementary diagrams" and triplet correlations have been treated as described 
in Ref. |T7. 
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The HNC equations also provide relationships between the correlation functions 
u n (ri, . . . , r n ) and the corresponding n-body densities. One of the quantities of primary 
interest is the pair distribution function g(r%, r 2 ) and the associated ( real-space ) static 
structure function 



S(r h r 2 ) = <y(n - r 2 ) + yj p 1 (r 1 )p 1 (r 2 )[g(r h r 2 ) - 1] . (2.5) 
The static structure function and the effective one-body Hamiltonian 

= -^TTT Vpi(r)v T^n (2 - 6) 

define the Feynman excitation spectrum through the generalized eigenvalue problem 

/^(n^Vi) = ^ / d 3 r 2 S(r 1; r 2 )^(r 2 ) , (2-7) 

which is readily identified with the inhomogeneous generalization^ of the well-known Feyn- 
man dispersion relational hcu(k) = h 2 k 2 /2m,BS(k). The states ip^\r), their associated 
energies hue, and the adjoint states 

1 H 1 (r)^\r) (2.8) 



are useful quantities for the impurity problem and for the representation of the dynamic 
structure function of the background film. 



B. The Static Impurity Atom 

The Hamiltonian of the N+l particle system consisting of N 4 He atoms and one impurity 



is 



ft 2 N 

+ UUro) + E ^(Iro - r,|) + H 



N 



(2.9) 



We adopt the convention that coordinate r refers to the impurity particle and coordinates 
Ti, with i = 1 . . . N to the background particles. Note that the substrate potentials E/ S ub(ri) 
and t/ s 7 ub (ro), as well as the interactions V^Qro — r 3 -|) and V(|rj — r 3 -|), can be different 
functions for different particle species. 

The generalization of the wave function fl2.2|) for an inhomogeneous iV-particle Bose 
system with a single impurity atom is 



l<i<j<N 



(2.10) 



exp 



Ki<N 



*iv(ri, ...,r N ) 



The energy necessary for (or gained by) adding one impurity atom into a system of N 
background atoms is the impurity chemical potential 
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// = Ejf +1 -E N . (2.11) 



Here, Ej^ +1 is to be understood as the energy expectation value of the Hamiltonian 
with respect to the wave function (|2.10|) . The further steps parallel those of the derivation 
of the background structure. 

The impurity density is calculated by minimizing the chemical potential ( p. 11] ) with 



respect to yp{(r ). This leads to an effective Hartree equation 

- ^ V o^ (r) (ro) + KL b (ro) + V H (v )] V ^(r ) = W P) (r ) (2-12) 

where Vff(ro) is an effective, self-consistent one-body potential for the single impurity. The 
lowest eigenvalue of Eq. ( |2.12j ) is the impurity chemical potential // = to, and the corre- 
sponding eigenfunction the density of the impurity ground state, J p[(r) = r/°)(r). 

In the systems studied below, translational invariance in the x — y plane is assumed, 
and the states are characterized by two quantum numbers, m and k m , associated with the 
motion perpendicular (m) and parallel to the symmetry plane (k m ). When unambiguous, 
as in the states r]^(ro) and </>^(ri), we shall use the single label (e.g. m) to collectively 
represent both quantum numbers. In particular, the states rf r ' (r ) depend only trivially on 
the parallel coordinate, 

^ r \r ) = ^ r \z )e ik ^ . (2.13) 

The unit volume is chosen as the size of the normalization volume. The corresponding 
energies are 

h 2 k 2 

t r = e r + —^, (2.14) 
2mj 



where e r are the eigenvalues of Eq. ( |2.12| ) for k\\ = 0. 



C. Impurity dynamics 



It is tempting to identify the higher-lying eigenstates of the "Hartree-equation" (|2.12 ) 



with the excited states of the impurity. This is legitimate only in a static approximation for 
the impurity features. However, such a simplification misses two important features: 

• If the momentum is a good quantum number, low-lying excited states can be discussed 
in terms of an effective mass. In our geometry, a "hydrodynamic effective mass" is 
associated with the motion of an impurity particle parallel to the surface; it is caused 
by the coupling of the impurity motion to the excitations of the background liquid. 
The local Hartree-equation ( p. 12 ) misses this effect. 



The effective Hartree-potential Vh{z) is real, i.e. all "excitations" defined by the 
local equation ( |2.12[ ) have an infinite lifetime. A more realistic theory should describe 



resonances and allow for their decay by the coupling to the low-lying background 
excitations of the host film. 
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Hence, a static equation of the type (|2.12j ) is appropriate for the impurity ground state 



only. The natural generalization of the variational approach to a dynamic situation is to al- 
low for time-dependent correlation functions u n (r , . . . , r n ; t). We write the time dependent 
variational wave function in the form 

<f>(t) = 1 e-^'Vtro, r 1; ...v N ; t) . (2.15) 

v I ^ ) 

Consistent with the general strategy of variational methods, we include the time dependence 
in the one-particle and two-particle impurity-background correlations, i.e. we write 

^(ro,!-!, ...r N ;t) = exp^ 5ui(r ;t) + ^2(1*0, rf, t) \E^ +1 (r , ri, r N ) . (2.16) 



Ki<N 



The time independent part remains the same as defined in Eq. (|2.10| ). The time-dependent 
correlations are determined by searching for a stationary state of the action integral 



S= [ C(t)dt 

d_ 



c(t)= miH 1 ^ - i n-jm , (2.i7) 



where H^ +1 is the Hamiltonian (|2.9| ) of the impurity-background system. 

The derivation of a set of useful equations of motion for the impurity have been given 
in Ref. EO. The final result is readily (and expectedly) identified with a Green's function 



expression, where the three-body vertex function describes an impurity atom scattering off 
a phonon, and is given in terms of quantities calculated in the ground-state theory. The 
motion of the impurity particle is determined by an effective Schrodinger equation of the 
form 



-,2 



2m 



ipi(r,uj) + J d 3 r' S(r, r', uo)ipi{r\ u>) = huipi{r, u>) , (2.18) 



where Vfl-(r) is the effective one-body potential of Eq. ( [2.12|) , and E(r, r', a>) is the impu- 
rity self-energy. Within the chosen level of the theory, S(r, t',uj) is describes three-body 
processes, 

where W mr (r) is the three-body vertex function that describes the coupling between an 
incoming 3 He particle to an outgoing 3 He in the state r as well as an outgoing phonon in 
state m. The detailed form of these matrix elements follows from the microscopic theory that 
has been described in length in Ref. |20|, it is not illuminating for the further considerations. 
The structure of Eqs. ( |2.18|) and ( [2.19| ) is of the expected form of an energy-dependent 



Hartree-equation with a self-energy correction involving the energy loss or gain of the im- 
purity particle by coupling to the excitations of the background system. It is the simplest 
form that contains the desired physical effects. 
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The energy denominator in Eq. ( 2.19|) contains the Feynman excitation energies defined 



in Eq. (|2.7| ) and the Hartree impurity energies of Eq. ( 2.12 ). These energies are too high, and 
we expect therefore that three-body effects are somewhat underestimated. A lowering of the 
spectra in the energy denominator by an impurity effective mass or by a more quantitative 
phonon/roton spectrum should have the effect of enhancing the importance of multi-particle 
scattering processes. Hence, it is expected that the binding energy of the surface resonance 
is still somewhat too high compared with experiments. On the other hand, it is not expected 
that a more quantitative spectrum in the self-energy should change the effective mass of the 
Andreev state considerably because the hydrodynamic backflow causing this effective mass 
is mostly caused by the coupling to ripplons, which are well described within the Feynman 
approximation. 

III. THE PHYSICAL MODELS 

We consider liquid 4 He adsorbed to a plane attractive substrate which is translationally 
invariant in the x — y plane, i.e. t/ su b(r) = U sn b(z). The systems under consideration are 
characterized by the substrate potential U 8n b(z) and the surface coverage 

n = dzpi(z) , (3.1) 
Jo 

where pi(r) = p\(z) is the density profile of the 4 He host system. This density profile is, along 
with the energetics, structure functions, and excitations of the film, obtained through the 
optimization of the ground-state (|2.2|) as outlined above; the procedure has been described 
in detail in Ref. [17. 



A. Ground state 

We have in this work studied the scattering properties of 4 He atoms for a number of se- 
lected substrate potentials and surface coverages; we have selected four cases for the purpose 
of a detailed discussion. The substrate potentials and the corresponding density profiles are 
shown in Figs. and [3[ The surface coverages are n = 0.3 A for each substrate potential; 
additionally we have considered the case n = 0.4 A for a Cs substrates as well as Mg for a 
case that is somewhat more attractive than the screened graphite, but also has a long range. 

Alkali metal substrate potentials are simple 3 — 9 potentials characterized by their range 
C3 and their well depth D. They have the form 



U suh (z) 



4Cf 
27D 2 



' ° 3 (3.2) 



z 9 z 3 



The range parameters C3 of these potentials have been calculated by Zaremba and Kohnii, 
the short-range z~ 9 term is phenomenological and fitted to reproduce the binding energies 
of a single atom on these substrates. Slightly more complicated is our model of a graphite 
substrate covered with two solid layers of 4 He. Most important for low-energy scattering 
properties is the coefficient C 3 of the long-range attraction, the values of C 3 for our substrates 
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of graphite, Cs, Na, and Mg are 180, 670, 1070 and 1750 K A 3 , respectively!!!!. The 
graphite- potential is relatively short-range but deep and produces a very visible layering 
structure of the background film; thus one obtains a rather "stiff" systemO. 

Fig. |2] provides a comparison of these four different potentials. It is seen that the alkali 
metal potentials are longer ranged, the magnesium substrate has the deepest potential well. 
At the opposite end of the potential strength is the Cs substrates. This substrate has 
received much attention in recent years because of the experimental finding that it is non- 
wettingiHiil. Note that the Cs-adsorbed films are metastable; they were examined with 
two purposes in mind. One is to generate a situation that is reasonably close the infinite 
half-space limit. Therefore, we have studied this case also for the larger surface coverage 
n = 0.4 A 2 . The second reason is that the nature of the low-lying excitations^ as well as 
that of the impurity states!! is somewhat different than those for the graphite model as will 
be seen below. 

The third case, a Na substrate, is an intermediate case which is of some interest for the 
nature of the 3 He bound states, whereas the Mg substrate is both deeper and longer ranged 
than the screened graphite. 



B. Background Excitations 

Our earlier workElSi! has discussed extensively the excitations of quantum liquid films 
adsorbed to various substrates. These studies have been concerned with the interpretation of 
neutron scattering experiments!! - !!, they have therefore focussed on excitations propagating 
parallel to the film. Typically, four types of modes were found: 

1. Surface excitations: At long wavelengths and on strong substrates, these are substrate 
potential driven modes with a linear dispersion relation 

uz{k) = c 3 k, (3.3) 

where C3 is the speed of third sound. At shorter wavelengths and in the case of 
an infinite half-space, the surface-mode is driven by the surface tension and has a 
dispersion relation 

u£(k) = —k 3 . (3.4) 

where a is the surface tension, and the density of the bulk liquid. In practice, 
the dispersion relation is linear only in a rather small momentum regime, and the 
ripplon dispersion relation ( |3.4[ ) is a quite good approximation!! for the surface — mode 

dispersion relation up to wavelengths of about 0.5 A . The theoretically predicted 
surface energy obtained from Eq. fl3-4|) by a k 3 ^ 2 fit to the dispersion relation is 

o —2 

o"th ~ 0.279 KA which compares favourably with the most recent experimental 
valueffl of a ex « 0.279 Kk~ 2 . 

2. Bulk Rotons: Films with a thickness of two or more liquid layers show already a 
quite clear phonon/roton spectrum. The spectrum starts at finite energy in the long — 
wavelength limit and contributes, in this momentum regime, very little to the strength. 
It takes over most of the strength in the regime of the roton minimum. 
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3. Layer Rotons: Films with a strongly layered structure also show excitations (identified 
as sound-like through their longitudinal current pattern) that propagate essentially 
within one atomic layer. These excitations have a two-dimensional roton with an 
energy below the bulk roton, and have been identified with a "shoulder" in the neutron 
scattering spectrum below the ordinary roton minimum. 

4. Interfacial Ripplons: On very weak substrates, like cesium, one can also have an 
"interface ripplonSH" . Its appearance can be understood easily from the following 
consideration: Consider first a film with two free surfaces. Obviously, this film would 
exhibit two ripplon modes, one at each surfaced. Now, a weak substrate is moved 
against one of the two surfaces. The character of the "ripplon" at this surface will not 
change abruptly; rather the circular motion of the particles will be somewhat inhibited, 
and the energy of the mode will rise. This is precisely what is seen in the energetics 
and the current pattern of this second mode on Cs. Stronger substrate potentials 
suppress this interface mode; to distinguish between an "interfacial ripplon" and a 
"layer phonon" one must look at the current pattern of the excitation!!. 

The above list of excitations is restricted to modes that can be characterized legitimately 
by a wave vector kn parallel to the surface. To calculate the response to particles impinging 
normally on the surface, one must also look at the types of excitations perpendicular to the 
surface. These cannot be rigorously classified by a wave number, but one should basically 
expect standing waves or resonances at discrete frequencies, approaching the excitations of 
a bulk system as the film becomes thicker. No ripplonic excitations or layer-modes should 
be visible in this case. 

The character of excitations is intelligently discussed by examining the dynamic structure 
function S(k, u). A general procedure has been developed in Refs. |T8| , |T§| to use time- 
dependent correlations for a quantitative calculation of the dynamic structure function. The 
simplest version of the theory is analogous to the Feynman approximationHS; the dynamic 
structure function in that approximation can be calculated directly from the solutions of 
Eq. 



S(k;u) 



(3.5) 



where the w (r) are adjoint states ( |2.8|) of the solutions of Eq. ( [2.7| ) for energy hu. The 
Feynman approximation has its well known deficiencies, and methods for its improvements 
have been derived which provide quantitative agreement with experiments. 

Previous work has concentrated on the theoretical interpretation of neutron-scattering 
experiments, it was therefore concerned with momenta parallel to the liquid surface. In the 
present situation we must allow for both parallel and perpendicular momentum transfer. 
We show in Figs. |] and |5] the dynamic structure function for parallel and p erpe ndicular 
momentum transfer. Fig. |] shows the picture familiar from previous workHiBEl: a low- 
lying excitation which can be identified with a ripplon by its dispersion relation and its 
particle motion, and a high density of states on the roton regime; note that the second 
lowest dispersion branch corresponds to the interfacial ripplon mentioned above. Note also 
that the modes below the continuum energy — /x 4 + ft 2 gn/2m4 are discrete; they have been 
broadened by a Lorentzian of the same strength to make them visible. 
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The situation is quite different for perpendicular scattering. Again, the discrete exci- 
tations below the the evaporation energy have been broadened. We see a dominant ridge 
basically along the dispersion relation of a Feynman phonon, and a high density of states 
in the regime of the roton. The ridge shows a number of "echoes" at shorter wavelengths; 
this is due to the finite-size of the film. But there are — expectedly — no excitations 
corresponding to the (interfacial) ripplons. 



C. Impurity Excitations 

Calculations of low-lying, bound states including the dynamic self-energy have been 



discussed extensively in Ref. |2(| we list here the most important ones demonstrating both 
the theoretical consistency as well as the quantitative reliability and highlight their relevance 
for scattering processes: 

• When applied to the bulk liquid, the ground state theory produces the correct chemical 
potentials of 3 He and hydrogenic impurities^. 

• In an inhomogeneous geometry, the static theory reproduces the binding energy of the 
Andreev state". The theory also predicts, even in its most primitive versionil, the 
existence of a surface resonance. 

• The dynamic theory predicts a hydrodynamic effective mass of the Andreev state of 
m* H /mj rs 1.35, to be compared to the value on 1.38 given by Higlev at a/.il somewhat 
larger than the value of m* H /mj « 1.26, reported by Valles et alM at the lower end 
of the value m* H /mj = 1.45 ±0.1 given by Edwards and SaamB In other words, our 
theoretical prediction is within the spread of experimental values. 

• The energy of the first excited surface state is lowered from about -2.2 K to -2.8 K, 
improving the agreement with the experimental valued of approximately -3.2 K no- 
tably. 

Similar to the obvious existence of interfacial ripplons, one also expects, on weak sub- 
strates, the appearance of an interfacial Andreev state. The binding energy of this state was 
found in Ref. ^ to be approximately -4.3 K. which is somewhat higher than the experimen- 
tal valued of -4.8 K. We attribute the difference to uncertainties in the substrate potential 
and the certainly oversimplified assumption of a perfectly flat surface. This state — being 
confined to a smaller area than the surface state — has always an energy that is higher 
than the Andreev state. Although it can in principle decay into a surface bound state, it 
has negligible overlap and hence its lifetime is practically infinite. With increasing potential 
strength, the energy of the substrate bound state increases; the state disappears completely 
on substrates somewhat more attractive than Na. Then, the "interfacial Andreev state" 
turns into a resonance to which a scattering particle can couple. Similar "resonances" can 
be found on Mg substrates even in the second layer; we shall return to this point further 
below. 

The two surface-bound states (and, if applicable, the interfacial Andreev state) can be 
described in the energy regime we are interested in reasonably well by an ti(k) = £j(0) + 
h 2 k 2 /2m^. Above the solvation energy of a 3 He atom, a sequence of impurity states can 
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exist that are spread out throughout the film; the detailed energetics of these states depends 
on the thickness of the film and the corrugation of the background liquid. 



IV. SCATTERING STATES 



The background and impurity excitations discussed in the previous section specify the 
possible energy loss channels for a scattering particle; we can now turn to the analysis of 
our results. 

The previous work has concentrated on the properties of bound impurity atoms, their 
effective masses, and the lifetime of resonances. Scattering processes are treated within the 
same theory, imposing asymptotic plane-wave boundary conditions on the solution of the 
effective Schrodinger equation ( [2.181 ): 

V>iO, r||) -> e ik ii' r [e~ ik±z + Re ik±z ] as z -f oo. (4.1) 

One of the key quantities of the theory is the elastic reflection coefficient R because 
it is directly influenced by the coupling of the motion of the impinging particle to the 
excitations of the quantum liquid. The absolute value of the reflection coefficient can differ 
from unity only if the self-energy S(r, r', uo) is non-hermit ian. This happens when the energy 
denominator in the self-energy ( |2.19| ) has zeroes; note that the quantum numbers m and 
r include both the motion of the particles parallel to the surface as well as the discrete or 
continuous degrees of freedom in the z-direction. 

Superficially, we appear to be describing a single-particle quantum mechanical scattering 
problem. In fact, a number of notions can be carried over from single particle models can be 
carried over, and simple phenomenological descriptions can be constructed at the level of a 
one-body theory. But the actual situation is far richer: Since the scattering film is composed 
of helium atoms, this is a generically non-local problem when viewed at the one-body level. 
Moreover, the film is dynamic: the incoming particle may produce excited states of the 
background. This may result in the capture of the particle and/or the emission of particles 
in states other than the elastic channel. 



A. Quantum reflection 



Generally, the amplitude of the wave function of an impinging particle of low energy is 
suppressed inside an attractive potential by the mismatch of the wave lengths inside and 
outside the potential if its range is small compared to the wavelength of the particle. As a 
consequence, the particles are almost totally reflected even if there is dissipation inside the 
potential (caused by the imaginary part of the self energy operator Q2.19| ) in our case) 



1 - \R\ oc k± as k± -> (4.2) 

and, consequently, s — > and ri ne i — > 0. The effect is called universal quantum reflectior^^. 

Quantum reflection can be described phenomenologically in an effective single particle 
picture with a complex optical potential. The many-body aspect of the problem is to de- 
termine the physical origin, the magnitude, and the shape as well as possible non-locality 
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of that "optical potential". The energy range where quantum reflection is visible in a 
many-body system like one of those considered here depends sensitively on the energy- 
loss mechanisms and calls for a quantitative calculation. Even in the limit of zero incident 
energy, the self energy ( |2.19 ) is non-hermitian, and thus allows in principle for sticking. Fur- 
thermore, this energy range is strongly affected by the long range features of the substrate 
potentialsilffi 

Specifically, in the 3-9 substrate potential models fl3T2]), the sticking coefficient s depends 
on the strength C3 of the potential: For a local potential with the asymptotic form C^z~ 3 as 
z — > 00, one can showS that the amplitude of the wave function inside the potential depends 
linearly on the normal momentum of the incoming particle. Increasing C3 makes the po- 
tential appear smoother for particles with long wave length, thus increasing the penetration 
depth and the probability to reach the film. Indeed, a calculation of the sticking coefficient 
from the non-Hermitian effective Schrodinger equation fl2.18|) gives, already in the distorted 
wave Born approximation (DWBA), s oc k. 

Inelastic scattering is, at low incident energies, only possible by coupling to ripplons. 
An analysis of the imaginary part of S(r, r', to) reveals that the contribution of the inelastic 
channels is proportional to E 7 ^ 2 which gives in the DWBA r; ne i oc E 4 . In other words, 
inelastic processes are negligible in the low-energy regime. 

Although it is not the main thrust of our paper, we have examined the low-energy reflec- 
tion probabilities. Fig. |6] shows three examples for the dependence of the sticking probability 
s ~ 1 — \R\ 2 on the the incident energy for normal incidence. While on graphite adsorbed 
films, quantum sticking is readily observable in the sense that the sticking coefficient starts 
to drop monotonically for wavelengths longer than 0.1 A , corresponding to energies less 
than 0.1 K, the linear dependence of s on k begins only at energies that are two to three 
orders of magnitude less for Cs adsorbed films (and similarly Mg and Na). 

Once the origin and properties of the optical potential for low-energy scattering are 
understood from a microscopic point of view, one may a posteriori construct simple, analytic 
models that provide, within the range suggested by the estimated accuracy of the microscopic 
picture, some flexibility to examine the dependence of s on features of the optical potential. 
A simple model consists of a local potential that approaches the substrate potential in the 
asymptotic region z — > 00 and that is approximated by a square well with a depth estimated 
from the binding energy of the Andreev state and a width of 15 A. The energy dissipation 
term can be included through a localized imaginary part of the typical magnitude of our self- 
energy. Such a model reproduces qualitatively the large values of s in the mK energy regime. 
Of course, the model fails to explain the dependence on ky , see Fig. [□]. For completeness, 
we should also add that retardation should be taken into account for quantitative results 
below 1-10 m K0. 



B. Ripplon coupling 

"Quantum reflection" as a generic phenomenon needs only some damping mechanism; we 
now turn to the task of many-body theory to identify and examine the physics that leads to 
damping. The basic physics is contained in the self-energy ( |2.19| ) used in our calculation; it 
includes the energy loss of an incoming particle with energy Tluj to a background excitation 
hu m , leaving the particle in the state t r . Within this model, damping is expected to be 
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somewhat underestimated because the possibility to emit two or more phonons has been 
neglected. 

Unless there is negligible overlap of the wave functions, the most efficient energy loss 
mechanism is the coupling to the lowest-lying excitation. These lowest lying excitations 
of the helium film are the surface waves (ripplons), hence one expects that the energy loss 
of the 3 He particle is dominated by the emission of a ripplon. This serves as a qualitative 
argument. However, the reality is more complicated for 3 He scattering because several states 
are accessible. The condition that an excitation contributes to the imaginary part of the self- 
energy is that the energy denominator of the self-energy fl2.19p vanishes, i.e. hu m + t r = hu, 
and there are several open channels even for vanishing incident energy. First, the particle 
can, although less efficiently, also couple to higher film excitations and can be promoted 
into either the second Andreev state or into a bound state in the bulk liquid. The reflection 
coefficients also depends visibly on the real part of the self-energy, and no quantitative 
statement can be made without proper treatment of both. The argument holds even at 
normal incidence, and infinitesimal asymptotic energy of the impinging particle. 

We show in Figs. |T|-|iO| a few typical examples of the reflection probability |-R(k||, k±)\ 2 
for scattering from 4 He films adsorbed on graphite, Na, and Cs substrates. In contrast to 
experiments on atomic scattering of 4 He from free 4 He surfaced^, there is evidently a strong 
dependence on the parallel wave vector kn which needs to be explained in terms the possible 
decay channels discussed above. Since it is unlikely that a specific feature is due to a delicate 
cooperation between film and impurity degrees of freedom, it is legitimate to discuss film- 
and single particle excitations independently. 

The fact that ripplon coupling is the dominant energy loss mechanism can be verified in 
various ways. The simplest one is the inspection of the self-energy (|2.19 ): The imaginary 



part of the self-energy is, with a few exceptions to be discussed below, localized in the 
surface region where the ripplon lives. The consequence is that, at energies below the roton, 
the wave functions of the impinging particle decays basically within the surface region. The 
effect can be seen in the wave functions and even better in the probability currents which 
basically decay within the surface region. A "resonance" in Figs. 12 and 13 will be discussed 
momentarily. 

From looking at Figs. [7|-|T0|, it appears that quantum reflection is seen only for the 
graphite substrate. As explained above, this is simply a consequence of the fact that the 
reflection becomes visible only at much lower energies on the alkali metal substrates. To 
demonstrate this, we have magnified in Fig. |TI| the low-energy region for the Cs substrate; 
consistent with Fig. || it is seen that the reflectivity starts to rise at impact energies of less 
than .001 K. 



C. Single particle resonances 

While the generic many-body aspect of all scattering and in particular damping mech- 
anisms must be kept in mind, one-body pictures can occasionally — as above for quantum 
reflection — provide useful paradigms in cases where the process under consideration can 
be described in terms of the degrees of freedom of a single particle. Such an effect is the 
coupling to single-particle resonances. A convenient and physically illustrative definition of 
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a resonance at an energy Hu is a large probability \ipi{r, u)\ 2 in the region of interaction. 
The resulting large dissipation will render |i?(a;)| 2 small. 

The peak of the wave function close to the substrate at k± ~ 0.4 A and z ~ 1.2 A 
shown in Fig. [12| is a very pronounced example of such a resonance. It displays exactly 
the phenomenon discussed above that the interfacial Andreev state turns into a resonance 
as the substrate strength is increased. The energy of this resonance is significantly reduced 
by the coupling to virtual phonons: The resonance has an energy of approximately 6 K in 
the static approximation ( |2.12| ). Including the dynamic self-energy corrections through the 
(real part of) £(co>), the resonance energy drops to approximately 1.3 K. The energy where 
the wave function has a strong peak in the vicinity of the substrate coincides with that of 
the dip in the reflection coefficient. Fig. |7] shows this for the special case of zero parallel 
momentum, but the agreement between the peak of the wave function and the minimum of 
the reflection probability persists at all parallel momenta. Also seen clearly in Fig. [1^ is the 
change in the phase of the wave as the resonance is crossed as a function of energy. 

The elliptic ridge of the reflection coefficient as a function of (k±, k\\) can be explained by 
the coupling of the interfacial Andreev resonance discussed above to the virtual excitations 
of the film. This has the consequence that the resonance acquires an effective massS rn* es . 
At zero parallel momentum, the position of the dip in the reflection coefficient agrees with 



the location of the resonance seen in Fig. 12. The shape of the ridge can be explained by 



assuming that all of the energy of the impinging particle is deposited in that resonance. 
Energy conservation and momentum conservation parallel to the substrate then leads to the 
relationship 



h 2 k\ h 2 k\ 

H - 

2m% 2m% 



1 _ JIjL 



(4.3) 



where e res is the energy of the resonance. Following the peak of the wave function in the 
resonance in the (k±, k») plane leads, within the accuracy that can be expected from such a 
relatively crude argument, to the same conclusion. Basically - and expectedly — the same 
resonances occur at other surface coverages; the precise location of the dip in the reflection 
varies due to the multitude of other open scattering channels. A similar resonance occurs 
in the more strongly attractive Mg substrates, the corresponding wave functions are shown 
in Fig. In this case, one finds a second resonance in the second layer which is, however, 
less pronounced. A list of energies and effective masses is given, together with the values for 
the Andreev state and the results of Ref. ^ of the bound states, in Table |. The effective 
masses were obtained by fitting the curve defined by Eq. (fO|) to reproduce the location of 
the peak of the wave function within the visible region. As pointed out above, the weaker 
substrate Cs has a bound state to which the scattering particle cannot couple, whereas the 
Na substrate is a marginal case. 

In all cases considered here we have found a significant dependence of the reflection 
coefficient \R\ on the parallel momentum, cf. Figs. |7-10|. Such a feature can not be explained 



within a local, complex single-particle model, is not also not seen in experiments on 4 He 
scattering off 4 He films/surfacesB. The feature is most pronounced on graphite and Mg 
substrates, cf. Fig. [F[ 

Also for the other substrates (see Figs. §, | and |10D, R(u) depends on the parallel com- 
ponent of the momentum. One sees similar, but broader ridges in the reflection coefficient, 
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but no sharp peaks in the wave function, cf. Fig. [14|. The effect can also be explained 
by the features of the impurity states within the film. But this time, the impurity states 
are not localized but are extended states that will, with increasing film thickness, develop 
to 3 He states dissolved in the 4 He liquid. Consistent with this picture, the energy of the 
resonances in Na and Cs adsorbed films decreases with increasing surface coverage n, until 
they become bound states, cf. [16|. This slipping below the threshold E = is best seen in 
the phase shift at E — > 0, which jumps whenever this occurs (in case of Cs: at n = 0.330 
and 0.380A" 2 ). 



D. Roton coupling 

The presence of roton excitations affects the scattering properties at two levels. First, 
at all energies, the coupling to virtual rotons is a significant contribution to the real "part of 
the self-energy; omitting these contributions by, for example, restricting the state sums in 
the self-energy (|2.19|) to energies below the roton minimum, leads to reflection coefficients 
that are, even at energies well below the roton minimum, about a factor of 2 smaller than 
when excitations in the roton regime are included. This is to some extent plausible since 
the roton is a reflection of the short-range structure of the system which is dominated by 
the core repulsion, and such effects should make the film look "stiffer" . 

At higher energies, the coupling to roton excitations also opens a new damping mech- 
anism. As seen in Fig. |^, "roton-like" excitations appear also for film excitations perpen- 
dicular to the surface, and the impinging particle can couple to these excitations. In our 
calculations, the effect of roton coupling should become visible at an energy of about 15 K, 
this is because we have used a Feynman-spectrum in the energy denominator of Eq. (|2.19 ). 
In previous work^l, we have scaled the energy denominators in the self-energy by an amount 
such that the roton is placed at roughly the right energy. We have refrained from this phe- 
nomenological modification since this procedure would also scale the ripplon away from its 
correct value which is already obtained in the Feynman approximation. 

Above 15 K, the film loses its elastic reflectivity for 3 He atoms completely. There is, of 
course, still the possibility of some inelastic scattering, but we consider this scenario unlikely 
from our experience with the propag ation of 'He impurities in bulk <He§. Hence, we expert 
that 3 He atoms will be completely absorbed by 4 He films when the impact energy is above 
the energy of the bulk roton. The effect is also seen quite clearly in the wave function of 
the scattering particle which does not penetrate into the film at all at energies above that 
of the roton. 



V. SUMMARY 

We have set in this paper the basic scenario for calculations of atom scattering processes 
from inhomogeneous 4 He. This work parallels similar research on scattering of 4 He atomsi! 
and also provides the groundwork for applications on the presently active area of atom 
scattering from 4 He clusters. 

Technically, the calculations presented here are somewhat simpler than those for 4 He 
atom scattering! 4 ! (since the Hartree impurity spectra appearing in the energy denomina- 
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tors in [2.19| are decoupled in parallel and perpendicular motion in contrast to the rip- 
plon/phonon/roton spectra) which enabled us to do a systematic study of the dependence 
of the reflection coefficient on the parallel momentum. 

When applicable, our general conclusions are very similar to those that we have drawn for 
4 He atom scattering: Most of the physics happens due to ripplon coupling, the wave function 
is substantially damped in the surface region. "Quantum reflection" does not come to bear 
until energies as low as 0.1 if on a graphite substrate, and .01 K or less on alkali metals. 
A second damping mechanism happening at higher energies is the coupling to rotons, this 
effect dampens the impurity motion completely; an equivalent effect is expected, and found, 
in bulk 4 He. 

A new aspect specific to 3 He scattering is the coupling to single-particle resonances within 
the film; such an effect will not be seen for 4 He scattering. We have demonstrated that the 
properties of the remaining reflected particles are directly influenced by the features of the 
impurity states within these films and that scattering experiments can directly measure the 
energy and the "effective mass" of these resonances. Fully acknowledging the experimental 
difficulty of the task, we hope that these findings will inspire further measurements on 3 He 
scattering off 4 He surfaces and films. 

Unfortunately it is difficult to make direct comparisons with experiments available 
todayEL One reason is that we are with our calculations apparently still too far from the 
bulk limit that a comparison is meaningful. This is most clearly seen in the oscillatory 
dependence of the reflection coefficient on the energy of the incoming particle even in a case 
if a relatively thick film without localized resonances within the film (Fig. |14|). There is 
also still pronounced non-monotonic dependence of the reflection coefficient on the surface 
coverage, cf. Fig. |16| . 

Further applications of our work are twofold: One is the application to scattering of 
hydrogen isotopes off 4 He surface. While experimental efforts in this area have been signifi- 
cantly strongerEJii^, the situation is less rich: The H impurity is only very weakly bound and 
can lose its energy only to the ripplon, in other words the imaginary part of the self-energy 
( P-19| ) comes from one state only. Moreover, the H atom must overcome a potential barrier 
of about 10 K to penetrate into the bulk liquid which makes the coupling to any interior 
degrees of freedom negligible. 

Similarly interesting is the possibility of scattering experiments of both 3 He and 4 He 
atoms off 4 He droplets. These experiments, to be carried out in the energy regime of a few 
tenth to a few degrees, would also couple to both surface and volume modes and should 
also more clearly display the coupling to excitations inside the droplets. In the spherical 
geometry, inelastic processes of the kind described here cannot occur at low energy because 
the continuous quantum number ky is replaced by the discrete angular momentum. Hence, 
all low-lying modes are discrete and the energy denominator of Eq. Q2.19 ) will normally 
be non-zero, in other words the self-energy is hermitean. Second, that these systems are 
not contaminated by substrate effects and should therefore allow for a cleaner interpretation 
of the results. We have learned that such scattering experiments have meanwhile been 
performed^ and show indeed the expected coupling of 3 He particles to roton-like excitations 
inside the droplets. 

Such experiments and calculations would also provide an ideal scenario to test ideas and 
procedures established in nuclear physics in a much better controlled and — in terms of the 
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underlying Hamiltonian — better understood physics. 

Calculations in this direction are in progress and will be published elsewhere. 
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TABLES 



TABLE I. Resonance energies and effective masses of the interfacial Andreev state on various 
substrates. The first line gives, for reference, the data of the Andreev state at the free surface, and 
the second the interfacial Andreev state on a Cs substrate (From Ref. |20| ). The last three lines 
give the results obtained here from scattering properties. Energies are given in K, the row labeled 
with "C" refers to the graphite plus two solid layers of 4 He model used in this work. 



Substrate 


Energy 


m*/m3 




-5.4 


1.3 


Cs 


-4.3 


1.7 


C 


1.3±0.3 


1.7±0.3 


Mg 


4.3±0.5 


1.6±0.2 
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FIGURES 






FIG. 1. The three classes of scattering channels are illustrated. The incoming particle can be 
(a) scattered elastically (top figure), (b) inelastically, (middle figure), or (c) adsorbed to the film 
(bottom figure). 
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FIG. 2. The figure shows the three substrate potentials for the films under consideration here: 
Graphite plus two solid helium layers (solid line), Mg (long dashed line), Na (short dashed line) 
and Cs (dotted line). 
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FIG. 3. The figure shows the four density profiles of the background liquid (solid lines) and 
the impurity location (long dashed lines) for which most of the present calculations were done. 
Graphite substrate results are marked with +-symbols, Na results with stars, and Cs results with 
crosses. Also shown is the interfacial Andreev state on a Cs substrate (short-dashed line marked 
with crosses). Coverages are n = 0.30 A for Cs, Na, and graphite, and n = 0.40 A for Cs. 
Profiles on Mg have been left out for clarity. 
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FIG. 4. Dynamic structure function S(k, uo) in Feynman approximation for a film with coverage 
of n = 0.400A"' on a Cs substrate and parallel momentum transfer. The solid curve shows the 
continuum boundary — /x + h 2 q?,/2m4 and the dashed line the bulk Feynman spectrum. 
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FIG. 5. Same as Fig. |] for momentum transfer perpendicular to the film. The horizontal solid 
line shows the continuum boundary — fx and the dashed line the bulk Feynman spectrum. 
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FIG. 6. Sticking on a graphite, Na, and Cs-adsorbed film of n = 300^4~ 2 . The square boxes 
in the upper left of the plot are the data of Ref. 0. Note that these data were taken at an impact 
angle of 60deg. 
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FIG. 8. Same as | for a Cs film of density n = 0.30A" 2 . 
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FIG. 9. Same as for a Na film of density n = 0.30A -2 . 
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FIG. 10. Same as g for a Cs film of density n = 0.40A" 2 . 
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FIG. 11. Fig. H is magnified into the regime of low k±_ to demonstrate that \R\ finally 
approaches unity. 



31 



N 




FIG. 12. The figure shows the wave function IVK^)! 2 of a 3 He as a function of distance z and 
perpendicular wave number p± for normal incidence. The left face shows, for reference, the density 
profile of the film and the back face the reflection coefficient R(p±)- The substrate is graphite plus 
two solid helium layers, the surface coverage is n = 0.300A . 
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FIG. 15. Reflection coefficient |i?| 2 for normal incidence on a film with n = 0.300A ' on a 
Cs-adsorbed film. The solid line is the result when all relevant intermediate states are kept in the 
state sums ( 2.19| ) whereas the dashed line is the result when the sum over intermediate states is 
truncated below the roton minimum. A new scattering channel opens at 3.5 K. 
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